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ABSTRACT 

We investigate the effect of possible a-priori "imprinting" effects of general 
relativity itself on satellite/spaceraft-based tests of it. We deal with some per- 
formed or proposed time-delay ranging experiments in the Sun's gravitational 
field. It turns out that the "imprint" of general relativity on the Astronomical 
Unit and the solar gravitational constant GM Q , not solved for in the so far per- 
formed spacecraft-based time-delay tests, induces an a-priori bias of the order of 
10~ 6 in typical solar system ranging experiments aimed to measuring the space 
curvature PPN parameter 7. It is too small by one order of magnitude to be of 
concern for the performed Cassini experiment, but it would affect future planned 
or proposed tests aiming to reach a 10~ 7 — 10~ 9 accuracy in determining 7. 

Subject headings: Experimental studies of gravity Ephemerides, almanacs, and 
calendars Lunar, planetary, and deep-space probes 
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PACS: 04.80.-y, 95.10.Km, 95.55.Pe 



1. Introduction 



Several space-based tests of the Einstein's General Theory of Relativity (GTR) 
have been performed, or attempted, in the more or less recent past by following an 
"opportunistic" approach, i.e. by suitably analyzing existing data sets of artificial satellites 
or interplanetary spacecrafts almost always built and launched for different original 
purposes (e.g. satellite geodesy, geodyna mics, planetolog y, etc.). A cornerstone result 



was the Cassini radio science experiment f lless et al.lll999l ) which lead to constraining the 
devia tion of the PPN parameter 7 from its gen eral relativistic value of unity at a 10~ 5 



level (IBertotti et al.ll2003l ; lAnderson et al.l 120041 ) through the measurement of the general 



relativistic time delay affecting the electromagnetic waves linking the Earth and the Cassini 
spacecraft during its journey to Saturn when they were in superior conjunction, i.e. aligned 
and on opposite sides of the Sun. Many other high-precision space-based tests, aiming to 



2008 



Ashbv et al. 


2009; 


Hobbs et al. 


2009; 


Turvshev 


2009 



Milani 



2009 ). Indeed, some 



theore tic al scenarios predict deviations from GTR at such a level (IDamour and Nordtvedt 
1993a|jbl; IDamour and Polyakovlll994l ; IDamour et al.ll2002l ); for a a recent overview see, e.g. 
Turvshev! fl2008h . 



In this paper we wish to critically discuss certain subtle issues pertaining the consistency 
of such data analyses interpreted as genuine tests of GTR. In Section [2] we outline the 
problem in general terms. We will move to concrete examples in Section [3] in which we will 
deal with the Cassini-like ranging experiments (Section 13. ip . Section H] is devoted to the 
conclusions. 



2. The issue of the a-priori "imprinting" of GTR itself in tests dedicated to it: 

general considerations 

In general, in such tests huge data sets from man-made interplanetary probes/satellites 
are confronted to given dynamical models of their motion in which the relativistic effect to 
be tested was explicitly included, with one or more solve-for parameters {P} accounting for 
it to be estimated in a least-square fashion along with many other ones {K}, not directly 
pertaining GTR. 

Of crucial importance for interpreting such data analyses as genuine tests of GTR is 
to clarify how the numerical values of the models' parameters {F} which have been kept 
fixed to certain reference values, i.e. which have not been solved-for, have been originally 
obtained. From the point of view of testing GTR, it is not enough that the resulting post-fit 



-4- 



residuals of a certain directly observable quantity are statistically compatible with zero at a 
good level; the standard data reduction procedure used for the original goals of the missions 
exploited for the GTR test considered may not be valid, in principle, for performing a truly 
unbiased, genuine check of GTR which is not a "tautology" . 

Indeed, if the primary task of a space-based mission is, for example, to reach a certain 
astronomical target with a given accuracy, the only thing that is important to this aim is 
that the dynamical models adopted to predict the probe's motion are accurate enough; 
this is usually quantitatively judged by inspecting the post-fit residuals of some directly 
measurable quantities like, e.g., ranges or range-rates. How the parameters {F} entering 
the models have been obtained, i.e. their a-priori values, does not matter at all: the only 
important thing is that the resulting fit of an existing set of observations is good enough to 
minimize the observable's residuals. 

Such an approach may, in principle, not be entirely adequate when the goal of the 
data analysis is testing a gravitational theory like GTR in an unambiguous, unbiased 
and self-consistent way. In this case, how the fixed parameters {F} of the models have 
been obtained does, in fact, matter. Indeed, if one or more of them {/} were previously 
obtained from different data of different bodies in such a way that they somehow retain a 
non-negligible a-priori "imprint" of the same effect we are now interested in, their use may 
bias the current test just towards the desired outcome yielding, for example, a very high 
accuracy confirmation. In this cases, it would be more correct to use, if possible, values 
of such "imprinted" parameters {/} which have been obtained independently of the effect 
itself whose existence we are just testing in the present data analysis, even if the accuracy of 
such different values of the "suspect" parameters {/} was worse. Alternatively, if, for some 
reasons, such "unbiased" values are not available, { /} should be included, if possible, in the 
list of the solved-for parameters along with the one(s) {P} accounting for the effect to be 
tested, and the resulting covariance matrix should be checked to inspect the correlations 
among them. The price to be paid may be an overall accuracy of the test not so high as 
that previously obtained^, but we would have more epistemo logically consistent, reliable 
and trustable tests. 



3. Application to some concrete cases 

3.1. The Cassini radio-science test and other proposed high-accuracy 

space-based measurements of 7 

To be more definite, let us look at the Cassini radio science test. In that case, the 
radiotechnical data of the spacecraft traveling to Saturn were contrasted with a set of 



In principle, certain confirmations may have their status changed. 
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dynamical models by JPL of its motion and electromagnetic waves propagation in such a 
way that a correction A7 to the GTR-p redicted value of the PPN parameter 7 was solved 
for, among other parameters, obtaining (IBertotti et al.ll2003l ) 



A7 = | 7 -1| = (2.1 ±2.3) x 10" 



other authors got ([Anderson et al.ll2004l ) 

A7 = I -y — 1 1 



-1.3 ±5.2) x 10" 



(2) 



Now, a physical parameter which is, of course, crucial in such a test is the gravitational 
constant GM & of the Sun, which is the source of the relativistic time delay put on the test. 
It was not estimated (IBertotti et al.ll2003l ; lAnderson et al.ll2004l ). so that its numerical value 
was kept fixed to the standard reference figure of the JPL DE ephemerides used. It does, in 
principle, contain an a-priori "imprinting" by GTR itself through the same effect itself that 
was just tested with Cassini, in particular by 7 itself. Indeed, the numerical value of GM & 
comes from the fixed value of the defining Gaussian constant^ 

k = 0.01720209895 au 3/2 d~\ (3) 

and from the value of the Astronomical Unifl not estimated in the Cassini tests, 

AU = 1.49597870691 x 10 11 (±3) m (4) 



through 



GM & = k' 2 AU 3 d" 2 = 1.32712440018 x 10 20 (±8 x 10 



m 3 s 2 . 



(5) 



AU was, in fact, obtained just through a combination of radar ranging of Mercury, Venus, 
and Mars, laser ranging of the Moon (making use of light reflectors left on the lunar surface 
by Apollo astronauts), and timing of signals returned from spacecraft as they orbit or make 
close passes of objects in the solar system (jStandishl 120041 ); thus, it is affected in a non 
negligible way, given the level of accuracy of the techniques adopted, by GTR itself and, 
in particular, by 7 which enters the PPN expressions for the time delay and bending of 
traveling electromagnetic waves. Thus, there exists, in principle, the possibility that the 
high-accuracy results of the Cassini radio science tests may retain an a-priori "imprint" of 
GTR itself through GM & (and the Astronomical Unit as well). 

Let us put our hypothesis on the test by making some concrete calculations; for the 
sake of clarity, we will refer to the Cassini radio science tests, but the conclusions may be 
considered valid also for any of the many proposed 7— dedicated missions. 



2 See on the WEB http://ssd.jpl.nasa.gov/7constants 



3 Here we will use au for the symbol of the Astronomical Unit, like m for the meter, while 
AU will denote its numerical value in m. 
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The GTR time delay experienced by electromagnetic waves propagating from point 1 
to point 2 is 

2R g , fn + ^ + r^ + R, 



At = — gin 11 ' z " g , (6) 

C V r l + r 2 - r 12 + itg/ 

where R g = 2GM Q /c 2 is the Sun's Schwarzschild radius in which G is the Newtonian 
gravitational constant, M Q is the solar mass, c is the speed of light in vacuum; n is the 
heliocentric coordinate distance to point 1, r 2 is the heliocentric coordinate distance to point 
2, and r 12 is the distance between the points 1 and 2. Eq. (ED is the expression actually used 
in the JPL'S Orbit Determination Program (ODP) used to analyze interplanetary ranging 
with planets and probes. In order to quantitatively evaluate the level of "imprinting" 
by GTR itself in the used value of the Astronomical Unit, let us assume r% equal to the 
Earth-Sun distance and let us vary r 2 within 0.38 au and 1.5 au to account for the ranging to 
inner planets; the maximum effect occurs at the superior conjunction, i.e. wheri@ ni ~ — n 2 , 
and ri2 ~ T\ + r 2 . It turns out that At ranging k4x 10~ 4 s, which is certainly not negligible 
with respect to the accuracy of the order of 10~ 8 s with which the light-time foi@ 1 au ta is 
actually measured ( jht t p : / /ssd . j pi. nasa . gov / ? const ants [ ) . As a consequence, the quantitative 
impact of the interplanetary ranging in the inner solar system to the determination of the 
Astronomical Unit is of the order 

dAU = cAt ranging = 1.14291 x 10 5 m, (7) 

not negligible with respect to the meter- level accuracy in measuring the Astronomical Unit; 
thus, gAU/AU= 8 x 1CT 7 . Differentiating eq. §5§ with respect to au and eq. ([7]) yield 

Thus, we conclude that the technique adopted to determine the numerical values of the 
Astronomical Unit and of the Sun's GM induced an a-priori "imprint" of GTR on them of 
8 x 10~ 7 and 2 x 10~ 6 , respectively. 

Let us apply this result to a typical radio science experiment in the solar system with 
r\ fixed to the Earth-Sun distance. By writing r 1 / 2 = Xi/ 2 au, with xx/ 2 expressing distances 



4 Here n = r jr. 

5 The value in km of the Astronomical Unit is obtained by measuring at a given epoch 
the distance between the Earth and a target body (a planet or a probe orbiting it) by 
multiplying c times the round trip travel time r of electromagnetic waves sent from the 
Earth and reflected back by the target body, and confronting it with the distance, expressed 
in AU, between the Earth an d the target bo dy at the same epoch as predicted by some 
accurate dynamical ephemeris fjStandishl 120041 ) . 



-7- 



in Astronomical Units, differentiation of eq. (|6]) with respect to au and GM Q , and eq. 
fl7])-eq. OH]) yield an "imprinting" effect of the order of 



5(At) 



At 



2 x 10" 



(9) 



GTR 



for r2 up to ten^l AU; it turns out that the largest contribution comes from dGM & . It is 
too small by one order of magnitude with respect to the performed Cassini radio science 
tests, but it should be taken into account in the future, more accurate experiments whose 
expected accuracy is of the order of 10~ 7 — 10 -9 , in the sense that the a-priori bias of GTR 
in the future determinations of deviations of 7 from unity will be as large as, or even larger 
than the effects one will to test, unless either GM & will be estimated as well along with 7 
itself or a value obtained independently of it will be adopted. 

About the first point, we mention that 7 and (3 were never been estimated in those 
solutions in which, among other things, also AU was fitted" instead, they were kept fixe d 
( Krasinskv and Brumberg 2004 ; Pitjeva 20051; Fienga et al. 2008 ; Pitjeva and Standish 



20091 : iFienga et al.l l2009h 7 Moreover, iFienga et all f l2009h used a modified version of 



their latest planetary ephemerides, named INPOP08d, in which they kept fixed AU and 
simultaneously fitted GM Q to the observations along with other parameters, but they did 
not simultaneously estimate the PPN parameters as well which were kept fixed. Instead, 
for the purposes outlined in this paper it would be important to simultaneously fit 7 and 
AU (or GM & ) and inspect the covariance matrix for the correlation among AU (or GM Q ) 
and 7. 

Concerning the last point, a possible choice may consist, in principle, of using a figure 
for GM & obtained from measurements of the solar gravitational redshift 



GMc 







RqC* 



(10) 



which, at first order, is independent of GTR0 itself. Latest measurements of the IR 
oxygen triplet 7772-7775, extrapolated to the Sun's limb ( jLopresto et al.lll99ll ). yield an 
accuracy of the orde r of a few percent level. Projects to imp rove it with future missions 
have been proposed (ICacciani et al.ll2006l ; iBerrilli et al.l 120091); f o r exam ple, the use of the 
Magneto- Optical Filter technique, developed by ICacciani et al.l ( 120061 ). would allow to 
reach a relative accuracy of 10~ 6 . Anyway, it must be noted that extracting GM & from the 
measured gravitational red-shift also requires the knowledge of the Sun's radius R Q , which 



6 The Cassini test was performed with r2 = 7 A3 au (IBertotti et al.ll2003l ). 

7 Indeed, it depends on the coefficient g 00 of the spacetime metric tensor; only the PPN 
parameter (3 enters g 00 with the term of order C(c~ 4 ), at present undetectable. 



-8 - 



is uncertain at a 10 4 level. Indeed, the commonly accepted value for the solar radius was 
for a long time f l A Ueru 1 1 9 73h 

#© = 695.99 Mm (11) 

(1 Mm = 10 6 m), although it was not clear how such a figure was obtained; moreover, 
no error bar was r eleas ed. Subseque nt observations of solar /—mode frequencies induced 
Schou et al. ( 19971 ) and Antia ( 1998 ). who used data from the Michelson Doppler Imager 
(MDI) mounted on the Solar and Heliospheric Observatory (SOHO) satellite and from the 
Global Oscillation Network Group (GONG) network, to conclude tha t the actual solar 
radius was 0.3 — 0.2 Mm smaller: indeed, the estimate by ISchou et al.l ( 119971 ) is 



R Q = 695.68 ± 0.02 Mm. 



In iBrown and Christensen-Dalsgaardl ( 119981 ) the value 

R Q = 695.509 ±0.026 Mm 



(12) 



(13) 



is reported from the Earth-based High Altitude Observatory's Sol ar Diameter Monitor 
campaign; it differs by about 0.5 Mm from the stand ard value by lAllenl (119731). being 
inconsistent with it at much more than 3 — a level. In iTakata and Goughl (j200l[ ) we find, 
from SOHO/MDI frequency data for p— mode frequencies, 



Rn = 695.69 ±0.14 Mm. 



(14) 



Note that the values by Brown and Christensen-Dalsgaard ( 1998 ) and Takata and Goughl 
( 2001 ) are mutually inconsistent, although at a 1 — a level only: the estimates by lSchou et al. 
( 1997 ) and Brown and Christensen-Dalsgaard _i 19981) are not consistent at 3 — a lev el. The 
/— mode and p— mode measurements (ISchou et al.lll997t ITakata and Goughl l200ll ) are, 
instead, consistent each other. By the way, it can be noted that the discrepancies among 
the different best estimates are larger than the associated errors; thus, as a conservative 
evaluation of the accuracy in knowing the solar radius we will assume such differences. 



4. Conclusions 

We have investigated the impact of possible a-priori "imprinting" effects of GTR itself 
on satellite/spacecraft data analyses specifically designed to test some general relativistic 
predictions. In particular, we considered the time-delay experiments conducted or proposed 
in the Sun's field with ranging to interplanetary spacecraft. 

Concerning the time-delay tests, the numerical values for the Astronomical Unit and 
the solar gravitational constant currently adopted in the ephemerides used so far to analyze 
spacecrafts' data retain an a-priori "imprint" of GTR itself of the order of 8 x 10~ 7 and 
2 x 10~ 6 , respectively. As a consequence, the bias in typical solar system radio-science 
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experiments is of the order of 10 -6 , which is one order of magnitude smaller than the 
accuracy level reached in the performed Cassini experiment, but it would be of concern for 
future planned tests aiming to measure deviations of the PPN parameter 7 from its general 
relativistic value at 10~ 7 — 10 -9 level. 

In order to have genuine, unambiguous and unbiased tests of GTR which are not 
"tautologic" , it would be necessary to either estimate the suspect parameters as well along 
with those accounting for the relativistic effect of interest or use values for them obtained 
independently from GTR itself. 
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